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Abstract 

If  gamma-ray  bursts  are  sources  of  ultra-high  energy  cosmic  rays, 
then  radiative  signatures  of  hadronic  acceleration  are  expected  in 
GRB  data.  Observations  with  the  Fermi  Gamma-ray  Space  Telescope 
(Fermi)  offer  the  best  means  to  search  for  evidence  of  UHECRs  in 
GRBs  through  electromagnetic  channels.  Various  issues  related  to 
UHECR  acceleration  in  GRBs  are  reviewed,  with  a  focus  on  the  ques¬ 
tion  of  energetics. 


1  Introduction 

Fermi  observations  of  GRBs  provide  a  new  probe  of  particle  acceleration  in 
the  relativistic  outflows  of  GRBs.  Some  generic  features  of  the  high-energy 
behavior  of  10  Large  Area  Telescope  (LAT)  GRBs  consisting  of  8  long-soft 
and  2  short-hard  GRBs  at  the  time  of  this  conference  have  been  identified, 
as  summarized  by  Ornodei  [T] .  These  are: 

1.  A  delayed  onset  of  the  ^  100  MeV  radiation  observed  with  the  LAT 
compared  with  the  start  of  the  keV /MeV  GBM  radiation; 

2.  Long-lived  LAT  emission  extending  well  after  the  Gamma  ray  Burst 
Monitor  (GBM)  radiation  has  fallen  below  background,  as  known  pre¬ 
viously  for  long  duration  GRBs  from  EGRET  [2j; 
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3.  Existence  of  a  distinct  hard  spectral  component  in  addition  to  a 
component  described  by  the  Band  function  in  both  long  and  short 
GRBs,  confirming  the  discovery  of  such  a  component  from  joint 
BATSE/EGRET  TASC  analyses  [3]. 

In  this  contribution,  the  possibility  that  these  features  can  be  explained 
by  UHECRs  in  GRBs  is  considered.  Because  protons  and  ions  are  weakly 
radiative  compared  to  electrons,  even  with  escaping  energies  E  ~  1020  eV 
needed  to  explain  UHECRs,  large  amounts  of  nonthermal  hadronic  energy 
are  required.  We  consider  energetics  arguments  to  constrain  models  for 
UHECRs  from  GRBs. 

2  Maximum  energy  release  in  GRBs 

Bohdan  Paczyriski,  in  his  seminal  article  on  hypernovae  [4],  proposed  an  ex¬ 
planation  for  the  large  apparent  isotropic  energies  of  GRBs  by  appealing  to 
the  enormous  energy  available  in  the  process  whereby  the  core  of  a  massive 
star  collapses  to  form  a  long  GRB.  He  noted  that  the  ~  10  M0  core  of  a 
Type  II  supernovae  carries  ~  5  x  1054  erg  of  rotational  energy,  which  could 
be  tapped  through  Blandford-Znajek  processes  if  surrounded  by  a  highly 
magnetized  torus  formed  during  the  stellar  collapse  event.  Paczynski’s  ar¬ 
guments  suggest  that  the  maximum  energy  available  from  core  collapse  su¬ 
pernovae  is  therefore 

£max  ~  1054  erg  ,  (1) 

with  perhaps  an  order-of-magnitude  less  in  the  coalescence  events  thought 
to  form  the  short  hard  GRBs. 

The  association  of  long-duration  GRBs  with  Type  Ib/c  rather  than  Type 
II  supernovae  does  not  alter  his  energetics  argument,  but  the  recognition 
that  GRBs  are  jetted  greatly  reduces  the  absolute  energy  requirements.  En¬ 
ergy  extracted  through  Blandford-Znajek  processes  is  likely  beamed  in  view 
of  the  rapid  rotation  of  the  newly  formed  black  hole.  Even  so,  we  can  still 
apply  the  energy  bound  given  by  eq.  ©  to  determine  the  plausibility  of 
models  with  large  energy  requirements. 

3  Energetics  of  UHECRs  from  GRBs 

The  idea  that  GRBs  could  accelerate  UHECRs  was  proposed  by  Waxman 
[5]  and  Vietri  0  in  1995,  and  subject  to  criticism  on  energetics  grounds 
though  here  related  to  whether  GRBs  within  the  GZK  radius  have  the 


requisite  volume-  and  time-averaged  luminosity  density.  We  can  reformulate 
their  argument  in  light  of  new  data.  The  Auger  results  show  that  the  differ¬ 
ential  energy  density  of  UHECRs  at  1020  eV  is  «  2  x  10-22  erg  cm-3.  The 
horizon  distance  for  1020  eV  UHECR  protons  is  ~  50  Mpc  (shorter  than 
the  mean-free-path  of  ~  140  Mpc,  because  this  is  the  distance  from  which 
protons  with  measured  energy  E  originally  had  energy  ~  2.7 E),  so  that  the 
required  emissivity  to  power  £  1020  eV  UHECRs  is  ecr  ~cx2x  10~22  erg 
cm~3/50  Mpc  «  4  x  1043  erg  Mpc-3  yr_1.  If  UHECRs  are  protons,  then 
ecr  ~  1044£44  erg  Mpc^3  yr_1,  with  £44  ~  1,  is  required  to  power  ^  1019 
eV  UHECRs,  noting  that  the  exponentially  decreasing  horizon  distance  with 
energy  approximately  cancels  the  oc  E~2  decrease  in  UHECR  energy  density 
over  this  energy  range. 

With  long-duration  GRBs  taking  place  about  twice  per  day  over  the  full 
sky  for  BATSE-type  detection  capabilities  [8j,  and  are  found  at  a  typical 
redshift  of  unity,  implying  a  local  space  density  for  unbeamed  sources  of 
~  2  x  365  yr_1C/[47r(4  Gpc)3/3]  «  0.3(C/0.1)  Gpc^3  yr_1,  consistent  with 
the  value  of  ~  0.5  Gpc-3  yr^1  found  in  more  detailed  treatments  [9j.  Here 
£  accounts  for  the  smaller  star-formation  activity  occurring  at  z  ~  1  than 
at  the  present  epoch  (lOj.  If  each  GRB  releases  an  amount  of  energy  Ecr  in 
UHECRs,  then  the  local  volume-  and  time-averaged  cosmic-ray  emissivity  is 
ecr  —  5  x  10_1°Tc’i?(C/0-l)  Mpc-3  yr^1.  Equating  this  with  the  emissivity 
required  to  power  the  UHECRs  implies  that  each  GRB  must  release  Ecr  ~ 
2  x  1053£44(C/0.1)  erg  in  UHECRs.  This  can  be  compared  with  the  typical 
electromagnetic  energy  release  per  GRB  of  4tt(4  Gpc)210~5J7_5  erg  cm~2  ss 
2  x  1052E-5  erg,  where  J-  =  10  5Jr_5  erg  cm  2  is  the  average  BATSE  long- 
duration  GRB  fluence.  Thus  the  energy  released  in  UHECRs  has  to  be 

10  x  the  energy  measured  from  electromagnetic  processes,  independent 
of  beaming. 

The  need  for  large  baryon  loads  in  GRB  blast  waves  is  confirmed  by 
detailed  fits  to  the  UHECR  energy  spectrum,  which  imply  a  factor  ~  10  -  100 
times  more  energy  in  UHECRs  than  observed  in  electromagnetic  radiation 
HD-  An  explanation  for  the  rapid  declines  in  Swift  X-ray  light  curves  in 
terms  of  escaping  UHECR  neutrons  from  photohadronic  production  also 
requires  highly  baryon-loaded  GRB  outflows  P2I- 

4  Minimum  bulk  Lorentz  factor 

A  ^-function  approximation  for  the  77  opacity  constraint,  given  the  detec¬ 
tion  of  a  7-ray  photon  with  energy  mec2e \  and  variability  time  tv,  implies  a 


minimum  bulk  Lorentz  factor 
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where  fe  is  the  vF„  flux  at  photon  energy  mec2e.  Thus  Tmin  ~ 
916[/_66i(3  GeV)/f„(s)]1//6,  where  the  uFl/  flux  is  lCT6/_6  erg  cm'2  s_1, 
using  values  corresponding  to  time  bin  b  for  GRB  080916C  [13].  This  GRB, 
at  redshift  z  =  4.35T0.15  and  luminosity  distance  cLl  =  1-25  x  1029  cm,  had 
a  total  10  keV  -  10  GeV  energy  fluence  F  =  2.4  x  10~4  erg  cm-2,  implying 
a  total  7-ray  energy  release  £ltiSO  —  8.8  x  1054  erg. 

Provided  the  target  photon  number  spectrum  has  an  index  softer  than 
—  1,  eq.  m  gives  a  result  within  ~  10%  of  a  numerical  integration  over 
spectral  parameters.  Issues  in  the  evaluation  of  the  uncertainty,  ATmjn,  in 
the  value  of  Tmin  include  (1)  uncertainties  in  redshift  and  vFu  spectral  flux; 
(2)  the  definition  of  variability  time  tv;  (3)  the  cospatial  assumption  that  the 
target  photons  are  made  in  the  same  region  as  the  high-energy  photon;  (4) 
the  assumed  geometry  and  dynamical  state  fH]  of  the  emission  region,  giving 
the  escape  probability  of  a  high-energy  photon;  (5)  statistical  fluctuations 
for  the  detection  of  a  high-energy  photon  that  furthermore  depend  on  the 
actual  high-energy  photon  spectrum. 

Writing  the  GRB  blast  wave  Lorentz  factor  T  =  qTmin,  a  consideration 
of  these  various  issues  show  that  q  ^  0.5  represents  a  reasonably  conservative 
expectation  for  the  actual  value  of  T.  The  energetics  of  a  given  model  can 
depend  strongly  on  T. 


5  Emission  radius:  internal  or  external  shocks? 

The  total  apparent  energy  release  of  GRB  080916C  is  £iSO  =  1055£,55  erg, 
with  £55  ~  1-  The  corresponding  deceleration  length  is  r^ec  =  1-2  x 
101'  /jrTg^3  cm,  where  103T3  is  the  coasting  Lorentz  factor  and  n(cm~3) 
is  the  external  medium  density.  The  implied  radius  for  internal  shock  emis¬ 
sion  is  r  =  r2ct„/(l  +  z)  =  6  x  1015r2L,(s)2  cm.  The  unexpectedly  large 
emission  radius,  close  to  the  deceleration  radius  when  q  ~  2,  has  led  a  num¬ 
ber  of  researchers  to  argue  that  the  LAT  radiation  originates  from  leptonic 
synchrotron  radiation  at  an  external  forward  shock  mm-  Another  pos¬ 
sibility  is  that  the  delayed  LAT  emission  results  from  upscattered  cocoon 
radiation  m • 


6  7  rays  from  UHECRs  in  GRB  blast  waves 


We  can  make  some  simple  estimates  to  deduce  the  total  energy  needed  to  ob¬ 
tain  bright  hadronic  7-ray  emission  from  GRBs  through  proton  synchrotron 
and  photopion  processes  pang  eq].  Here  we  consider  only  UHECR  proton 
acceleration,  using  parameters  appropriate  to  GRB  080916C. 

The  internal  photon  energy  density  ti7  =  d\<&  /  R2T2c,  where  = 
1(T5$_5  erg  cm-2  s^1  is  the  measured  energy  flux,  R  =  T2ct/(1  +  z)  is  the 
shock  radius,  and  t  is  a  fiducial  timescale  (corresponding  to  tv  for  internal 
shocks,  or  the  GRB  duration  for  an  external  shock).  Writing  the  magnetic- 
field  energy  density  u'B  =  esu'  where  u'B  =  B2/8it1  then  the  magnetic 
field  in  the  emission  region  of  GRB  080916C  is  B(kG)  =  2.0\/eB<&-5 /rp(s) • 
The  Hillas  criterion  whereby  the  Larrnor  radius  rf  =  mpc2n('/eB  <  A R'  = 
R/T,  where  A R'  is  the  comoving  shell  width,  implies  that  the  escap¬ 
ing  UHECR  proton  Lorentz  factor  =  T^'p  7  ceBT2t/(  1  +  z)mvc2.  For 
GRB  080916C,  this  relation  implies  that  UHECRs  can  be  accelerated  to 
7 p  7,  (2dLe/Tmpc2)y/2/KeB^ / c  =  4  x  lO^y/eB&s/Ts,  independent  of  time 
(note  that  B  oc  t-1,  with  the  divergence  at  t  — >•  0  prevented  due  to  the  time 
dependence  of  <f>),  or  to  energy  Ep  7  2  x  1021R(kG)Tp(s)  eV. 

A  further  restriction  on  the  maximum  proton  energy  is  obtained  by  bal¬ 
ancing  the  acceleration  rate,  given  by  c/(</>rf),  with  <p  /§>  1,  with  the  syn¬ 
chrotron  loss  rate.  This  gives  7 sat,P  =  2  x  lO^Ts/p (^/lO)S(kG),  compa¬ 
rable  to  the  value  obtained  above  from  the  Hillas  condition  for  the  chosen 
parameters. 


6.1  Proton  synchrotron  energy  requirements 

The  proton  synchrotron  energy  loss  timescale,  as  measured  by  an  observer, 
is  tsyn  =  3mec2(l  +  z)/[4r/i3ccJTW^7p],  where  fi  =  rne/rnp.  The  typical 
photon  energy  (in  mec2  units)  of  the  measured  proton  synchrotron  emission 
is  £p,Syn  —  r/iR7;2/[(l  +  z)Bcr\,  where  Bcr  =  4.414  x  1013  G  is  the  critical 
magnetic  field.  From  this,  one  obtains  the  jet  power  associated  with  the 
magnetic  Held,  given  by 


Lb  = 
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m-  The  absolute  energy  requirements  are  £ai }S  —  LBtsynfb/(l  +  z)  = 
1.6  x  1059T36//3/f,  t%n(  10  s)/E7(100  MeV)2/3  erg,  where  j\,  is  a  beaming  fac¬ 
tor.  In  the  scenario  of  Ref.  [IS],  the  delayed  onset  corresponds  to  the  time 


for  protons  to  accumulate  and  cool  such  that  they  are  radiating  most  of  the 
proton-synchrotron  photons  near  energy  E 7 

Such  large  energies  disfavor  this  interpretation  for  the  LAT  radia¬ 
tion.  To  salvage  the  proton  synchrotron  model  [IB],  a  narrow  jet  open¬ 
ing  angle  of  order  1°  along  with  a  value  of  q  =  0.5  gives  £abs  = 
4  x  1053(r3/0.5)16/3(/b/10~4)  t%n(  10  s)/£'(100  MeV)2/3  erg,  within  accept¬ 
able  ranges.  Another  way  to  avoid  the  large  energy  losses  is  to  suppose  that 
the  UHECR  protons  are  accelerated  to  their  allowed  maximum  energy  and 
radiate  proton  synchrotron  photons  that  cascade  to  energies  ^  100  MeV. 
This  possibility  seemed  unlikely  given  that  the  spectrum  of  GRB  080916C  is 
consistent  with  a  single  Band  function  m-  The  detection  of  distinct  com¬ 
ponents  in  GRB  spectra  suggests  that  a  cascading  interpretation  be  more 
carefully  considered  mm- 


6.2  Photohadronic  energy  requirements 

Proton-photon  interactions  making  secondary  pions,  7  rays,  and  neu¬ 
trinos  represents  another  likely  channel  for  making  a  7-ray  component 
identifiable  in  the  Fermi  data.  The  efficiency  for  extracting  the  energy 
of  a  proton  with  escaping  energy  Ep  from  photohadronic  processes  can 
be  written  as  r]pi(Ep )  =  tdyn/tpl(Ep )  =  (R/Fc)tpJ(Ep),  where  tdyn  is 
the  dynamical  time  scale,  the  rate  for  photohadronic  energy  losses  is 
tz}(Ep)  =  c{Kin(Jrn)  [7°  de'n'(e '),  and  the  comoving  photon  spectrum 

n'{e')  =  d2Lfe/(mec3e'  2R2T2),  with  e'  =  (1  +  z)e/T.  Here  Kpiapl  =  70/rb 
above  threshold  photon  energy  defined  by  7V  =  400  HU.  Defining  protons 
with  energy  E^k  which  interact  at  threshold  with  photons  with  energy  epk 
at  the  peak  f€pk  of  the  vFu  spectrum,  then  E^k  =  400mpc2T2/[(l  +  z)epk  = 
7  x  1016r2/epfc  eV.  One  obtains 


Kpiapid2Lftpk 
r4mec4f„(l  -  b)epk 


/— 6 
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[5] [221  [23].  Here  b(<  0)  is  the  vFu  index  above  epk,  so  that  fe  =  ftpk{e/ePk)b 
when  e  >  epk,  and  rjpy(Ep )  =  r)p~f(Elpk)(Ep/ E'^)1^ .  Likewise  in  the  asymp¬ 
totic  limit  e  <C  epk,  r]pi(Ep )  =  r/pi(EPk)(Ep/ E^)1-0, ,  where  a  is  the  uFu 
index  at  energies  below  the  peak  energy  (provided  a  >  1). 

The  total  efficiency  for  photohadronic  production  depends  on  the  spec¬ 
trum  of  accelerated  protons  and  ions.  If  an  E p2  spectrum  is  assumed,  then 
the  efficiency  is  reduced  in  proportion  to  the  number  of  decades  of  weakly 
radiating  low-energy  protons.  This  might  not  be  too  severe  if  the  lowest 


energy  protons  have  escaping  energy  Ep  ~  T2mpc2  ~  Tg  PeV.  Furthermore, 
nonlinear  shock  acceleration  in  colliding  shells,  and  second-order  processes 
in  the  shocked  fluid  can  give  harder  cosmic-ray  spectra,  though  consistency 
when  fitting  the  measured  UHECR  spectrum  would  constrain  the  assumed 
accelerated  cosmic-ray  spectrum. 

The  final  expression  in  eq.  (HJ),  specific  to  parameters  of  GRB  080916C, 
shows  that  ~  1  -  10%  of  the  energy  can  be  extracted  via  photohadronic  pro¬ 
cesses  by  protons  with  energy  Ep  ~  E^k ,  with  more  than  half  of  this  energy 
released  in  the  form  of  leptons  and  photons  which  generates  an  electro¬ 
magnetic  cascade  emerging  in  the  form  of  7  rays  when  the  system  becomes 
optically  thin  to  77  processes.  At  Ep  E^k,  a  larger  fraction  is  extracted. 
Consequently,  the  energy  of  baryons  must  be  ~  10  to  100  x  greater  than  the 
nonthermal  lepton  content  if  a  comparable  amount  of  electromagnetic  radia¬ 
tion  is  to  be  emitted  from  hadronic  as  leptonic  processes.  Note  the  sensitive 
T-dependence  of  the  photohadronic  efficiency,  i^pl{E^k)  oc  T~4,  by  compari¬ 
son  with  r77(e)  25  T"6.  The  implications  for  neutrino  and  7-ray  production 
are  considered  in  [2Tj. 

7  Summary 

In  this  contribution,  we  have  sketched  the  energy  requirements  for  GRBs  to 
be  sources  of  UHECRs,  and  for  electromagnetic  signatures  of  ultrarelativis- 
tic  hadrons  to  be  found  in  the  Fermi  data.  The  large  amounts  of  energy 
needed  has  been  noted  many  times  in  the  past,  whether  from  proton  syn¬ 
chrotron  [25],  26]  or  photohadronic  processes.  Although  the  large  energy 
requirements  make  uncomfortable  demands  on  7-ray  emission  models,  an 
internal  consistency  is  found  insofar  as  the  baryon  load  in  long  GRBs  must 
be  large  given  their  relative  rarity  within  the  GZK  radius,  and  that  enor¬ 
mous  energies  are  available  from  the  rotational  and  accretion  energy  in  the 
newly  forming  black  holes.  Future  Fermi  observations  and  the  possibility  of 
detecting  PeV  neutrinos  from  GRBs  with  IceCube  could  establish  whether 
GRBs  are  the  sources  of  UHECRs. 
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